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ABSTRACT 

Synchrotron emission is believed to be a major radiation mechanism during gamma-ray bursts' 
(GRBs) prompt emission phase. A significant drawback of this assumption is that the theoretical 
predicted spectrum, calculated within the framework of the "internal shocks" scenario using the 
standard assumption that the magnetic field maintains a steady value throughout the shocked region, 
leads to a slope F v oc v^ 1 / 2 below 100 keV, which is in contradiction to the much harder spectra 
observed. This is due to the electron cooling time being much shorter than the dynamical time. In 
order to overcome this problem, we propose here that the magnetic field created by the internal shocks 
decays on a length scale much shorter than the comoving width of the plasma. We show that under 
this assumption synchrotron radiation can reproduce the observed prompt emission spectra of the 
majority of the bursts. We calculate the required decay length of the magnetic field, and find it to be 
~ 10 4 — 10 5 cm (equivalent to 10 5 — 10 6 skin depths), much shorter than the characteristic comoving 
width of the plasma, ~ 3 x 10 9 cm. We implement our model to the case of GRB050820A, where 
a break at < 4keV was observed, and show that this break can be explained by synchrotron self 
absorption. We discuss the consequences of the small scale magnetic field scenario on current models 
of magnetic field generation in shock waves. 



Subject headings: gamma rays: bursts — gamma rays: 
radiation mechanisms: non-thermal 



theory — magnetic fields — plasmas 



1. INTRODUCTION 

A widely accepted interpretation of the non-thermal 
radiation observed during the prompt emission phase 
of gamma-ray bursts (GRBs) is that synchrotron emis- 
sion is a leading radiation mechanism during this 

1993iJbt 



phase (|M eszaros et ql"1 IT993l IMeszaros fc R,ees I 

iKatzl 119941 IRees fc Meszaros I 



Meszaros et al I 11994 L, 

11994 iTavani Ill996aj) . Indeed, early works found that 
the majority of bursts show spectral slopes in the 
~ 1 - 200 keV r a nge of vF„ oc v a , with a ~ 4/3 
llTavani I Il996albt iCohen et al. I Il997t ISchaefer etaT\ 
1998; Fro ntera et al. II200C1) . which is in accordance with 
the predictions of the optically thin synchrotron emission 
model, provided that the synchrotron cooling time of the 
radiating electrons is longer than the emission time. In 
addition, recen t comprehensive analysis of the brightest 
BATSE bursts l|Preece et al. ll2000HKaneko et al. Il200l 
found that the distribution of the low energy spectral 
slope peaks at a ~ 1, and a significant fraction of the 
bursts show spectral slope consistent with a ~ 4/3. 

The idea that synchrotron emission is the leading 
radiation mechanism had gained further support by 
modeling the more detailed observations of the after- 
glow phase in GRBs, which are foun d to be in good 
agreement with this mo del prediction (iSari et al. Ill996t 



11992111994 ISari fc Piranlll997ft . magnetic fields are gen- 
erated by shock waves. Electrons are accelerated to high 
energies by the same shock waves, which thus provide 
the necessary conditions for synchrotron radiation. The 
mechanisms of energy transfer to the magnetic field and 
to accelerated electrons are not fully understood. It is 
therefore common to parametrize the energy densities 
in the magnetic field and in the energetic electrons as 
fractions cb and e e of the post shock thermal energy, 
where the values of e e and £b are inferred from ob- 
servations. By modeling GRB afterglow emission data 
(TWiiers fc Galama 1 119991: iFreedman fc Waxman 1 120011: 
iPanaitescu fc Kumar 1120011 12002|) . the parameters val- 
ues are found to be at least few percents in most 
of the cases, and in some cases close to equipartition 
(TWiiers fc Galama IH991 iFrail et al. ll200?ilh 

The fact that a significant fraction of the bursts 
show spectra that are too hard to account for i n the 
optically thin sync hrotron model llCrider et al. I Jl997t 
iPreece et all Il99l 12001 iGhirlanda et al. 1 12003|) had 
motivated works on alternative emission models. These 
include Synchrotr on self Comp t on (SSC) sca t tering , 
first suggested by iLiang I lll997fl: ILiang et al. I l)1997|) . 
Compton drag l|Lazzati et al. I 120001) . upscattering of 



Meszaros & Rees 1997; Waxman 1997a b; Sari et al. 
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Panaitcs cu fc Meszaros 1 11998: Wii ers fc Ga lama 



synch r otron self absorbed photo ns llGhiscllini fc Celotti 
1999; IPanaitescu fc Meszaros! 120001: iKumar et al. 



In the standard internal/external shock scenario of 
GRBs (the "fireball" model scenario; IRees fc Meszaros! 
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20061) , and Compton sc atterin g of photospheri c pho- 
tons llMeszaros fc ReeTI 120001 IMeszaros et al. 1 120021: 
iPe'er et al. 1120051 1200(1 and references therein). While 
in principle these models can reproduce a hard spectral 
slope, common requirements to all models involved 
inverse Compton scattering as a leading radiation 
mechanism is that at the emission radius the optical 
depth to scattering is high, and that es/e e ^ 1. This 
last requirement, in turn, can lead to extensive radiation 
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at very high MeV) energies, thus to low radiative 
efficie ncy at the sub-MeV energy range ijDerishev et al. I 
1200 1|) . Additional drawback of SSC models is the wider 
spread in the peak energy distribution (compared to 
synchrotron model results for similar range of param- 
eter dispersion) , which might not be co nsistent with 
the data (e.g., iZhang fc Meszaros 2002J). Therefore, 
these models put various constraints on the allowed 
parameter space region during the emission phase (e.g., 
IZhang fe Meszaros II200H iPe'er fe Waxman Il2004|) . 

An argument raised a gainst the synchrotron mech- 
anism | G hisellini et al. I l|2000D : see also discussion in 
IZhang fe Meszaros I i)2004T) ] is that the inferred values of 
the free model parameters, in particular the strength of 
the comoving magnetic field during the prompt emission 
phase, B' ~ 10 5 — 10 6 G, implies that the radiating elec- 
trons are synchrotron cooled much faster than the dy- 
namical time. This, in turn, leads to a spectrum with 
slope vF v oc v 1 / 2 below ~ 100 keV, which is in conflict 
to the much harder spectra observed in this energy range. 
In order to overcome this problem, it was suggested 
that the energy distribution of radiating electrons has 
a smooth cutoff, and that the pitch a ngles of these elec- 
trons are anisotropicallv distributed l|Llovd fc Petrosianl 
1200(11 ILlovd-Ronning fc Petrosian 112002ft . 

A crucial underlying assumption in this analysis, is 
that electrons radiate on a length scale comparable to the 
entire comoving width of the shocked plasma. For plau- 
sible assumptions about the number density and char- 
acteristic Lorentz factors in GRBs, this assumption can 
only hold if the magnetic field maintains approximately 
const ant value on a scale of w 10 9 skin depth ijPiran I 
120051) . 

Generation of magnetic field in strong, relativistic 
shock waves is still poorly understood. Two stream insta- 
bility of flow past shoc k waves can, in principle g enerate 
strong magnetic field ijMedvedev fc Lo eb ||199 9[). How- 
ever, state of the art numerical models iSih/neFal. I2003t 
iFrederiksen et al. 120041 iNishikawa et al. I2005|) can only 
trace the evolution of this field on a characteristic scale 
of few tens of skin depths at most, due to the huge 
numerical effort involved. The evolution of the mag- 
netic field on larger scale therefore still remains an open 
question. While some models predict that the magnetic 
field saturates at a v alue close to equipartition (e.g., 
Uaroschek et al. P2004T). several authors find m uch weaker 
magnetic field ijWiersma fc Achterberg 120 04) . or argued 
that the creat ed magnetic field quickly decays by phase- 
space mixing ijGruzinov 1120011) . 

Motivated by these uncertainties o n the length scale 
of the magnetic field, iRossi fc Rees I l|2003j) suggested a 
model for GRB afterglow emission, in which the magnetic 
field decays on a length scale shorter than the shocked 
region scale. In that work, however, the decay length of 
the magnetic field was not specified. 

In this paper we show that by assuming that the mag- 
netic field decays on a length scale shorter than the co- 
moving scale, the observed prompt emission spectra of 
many GRBs can be reproduced. Thereby, this assump- 
tion allows to overcome the "fast c ooling time" problem 
inferred bvlGhiscllini et al. (2000]). We calculate in sec- 
tion 33 the values of the free model parameters that can 
account for the GRBs prompt emission spectra. We show 



that the decay length of the magnetic field that is consis- 
tent with the observed spectra is ~ 10 4 — 10 5 cm, which 
is i=s lO 5 5 skin depths. We then implement our model in 
section <J2]to the specific case of GRB050820A, where a 
low energy break at ~ 4 keV was observed. We summa- 
rize our results and discuss the implications of our model 
in view of current models of magnetic field generation in 
relativistic shock waves in section 21 

2. THEORY OF SMALL MAGNETIC FIELD LENGTH 
SCALE: CONSTRAINTS ON MODEL PARAMETERS SET 
BY OBSERVATIONS 

We adopt the framework of the internal shock scenario 
and assume that variability in the Lorentz factor T of the 
relativistic wind emitted by the GRB progenitor leads to 
the formation of shock waves within the expanding wind 
at r adii much larger than the u nderlying source size (see, 
e.g.. IZhang fc Meszaros 1 12002(1 . We assume that these 
shock waves, produced at characteristic radius r from 
the progenitor, are the source of the magnetic field. We 
introduce a new length scale Ar' B , which is the comov- 
ing length scale characterizing the decay of the magnetic 
field. This decay length is much shorter than the co- 
moving width of the plasma Ar' ~ r/T. We derive in 
this section the constraints on the model parameters as 
inferred from observations. 

The radiating electrons are accelerated by the shock 
waves to a power law distribution with power law index p 
above some characteristic energy 7 m i n m e c 2 . Synchrotron 
radiation by these electrons is the main emission mecha- 
nism, therefore the break energy observed in many bursts 
at ej™- > 100 keV is attributed to synchrotron radiation 
from electrons at 7 m i n . Denoting by j c the Lorentz factor 
of electrons that cool on a time scale equal to the dynam- 
ical timescale, and by e° b ' the characteristic observed en- 
ergy of photons emitted by synchrotron radiation from 
these electrons, the requirement that the spectral slope 
vF v oc v a has a characteristic spectral index a ~ 4/3 be- 
low - 100 keV leads to ef- > 100 keV. The value of ef- 
can not be much greater than 100 keV, in order to ensure 
high radiative efficiency (see discussion in 21 below). 

The requirement that the spectral slope is not harder 
than 4/3, as is the case in a significant fraction of the 
bursts, implies that in these bursts inverse Compton scat- 
tering and thermal emission component do not play a 
significant role in producing the spectra below 100 keV. 
These conditions can be translated to the demand that 
the emission radius r is larger than the photospheric ra- 
dius, r p h. Additional two constraints are that the ob- 
served flux vF v and the synchrotron self absorption en- 
ergy ef sa , which produces a low energy break, are con- 
sistent with observations. 

The observational constraints can therefore be written 
as a set of equations in the form: 

(a) e% > 100 keV, 

(6) ef- > 100 keV, 

(c) r >r ph , (1) 

(<£) vF° h - ~ 10" 7 ergs" 1 cm" 2 , 

(e) ef sa <1 keV, 

We now apply the set of equations JIJ describing the 
constraints set by observations to constraints on the un- 
certain values of the free model parameters. 

Assuming variability in the Lorentz factor AT/T ~ 1 
on timescale At of the expanding relativistic wind, shocks 
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develop at radius r ~ T 2 cAt. Due to Lorentz contrac- 
tion, the comoving width of a plasma shell is Ar' ~ r/T. 
We use the standard fireball model assumption, in which 
the burst explosion energy is initially converted to kinetic 
energy. For isotropically equivalent central engine lumi- 
nosity L which is time independent over a period At, the 
isotropically equivalent number of protons ejected from 
the progenitor during this period is N p w LAt/Tm p c 2 . 
Therefore, the comoving number density of protons in 
the shock heated plasma is given by 



4:TTr 2 T 2 c 



= 1.8 x 10 13 L r^ 2 T^ 2 Co cm" 



n' p {r) 

(2) 

where £ is the compression ratio (£ — 7 for strong 
shocks) and the convention Q = 1Q X Q X is adopted in 
CGS units. Assuming that the proton internal energy 
(associated with the random motion) in the shocked 
plasma is 9 p m p c 2 , the comoving internal energy density 
is u' = n p p m p c 2 . The value of 6 p is not expected to be 
much larger than a few at most for mildly relativistic (in 
the comoving frame) shock waves. The magnetic field 
carries a fraction es of the internal energy density, thus 
the comoving magnetic field strength is given by 



B' = ^8tt€ B u' 

= 4.6 x io 5 l 1/2 r^ 1 r; 



-11/2 



B 



1/2 ,1/2 



,-0.5 °p,0 SO 



G. 



(3) 



We assume that a fraction e p i < 1 of the electron pop- 
ulation is accelerated by the shock waves to a power law 
energy distribution with power law index p above 7 m i n 
(and below 7 max )- Assuming that a fraction e e of the 
post-shock thermal energy is carried by these electrons, 
the minimum Lorentz factor of the energetic electrons is 
given by 



7n 



€pl 



V m ° ) lo. 



(7min/7max) 



X \&(p) 



-0.5 e„ 



:,o u p 



,o*(p), 



(4) 



where characteristic value log(7 max /7 m i n ) — 7 was used. 
The function ^(p) determines the dependence of the 
value of 7 m in on the power law index p of the acceler- 
ated electrons, and is normalized to ^(p = 2) = 1. A 
full calculation of this function for various values of the 
power law index p is given in appcndix lA.il Using equa- 
tions |3 and the break in the spectrum from burst at 
rcdshift z is observed at 



^ob. 



3 5/2 ,l/2, Tf2 , 



e - 2 e 1/2 (5) 
-0.5 fc P ;,o fc B,~0.5 v ' 



x^oCo /Z * 2 WkeV. 

Electrons in the shocked region propagate at veloc- 
ity close to the speed of light. Therefore, electrons 
cross the magnetized area in a comoving time k, Ar' B / c. 
Since this is the available time for electrons to radiate, 
equating the synchrotron cooling time and the crossing 
time of this area gives the cooling break of the elec- 
trons energy distribution, which occurs at Lorentz factor 
7 C = (9m 3 c 6 )/(4g 4 B' 2 Ar^). Photons emitted by elec- 
trons at 7c are observed at energy 

97 

B,7 u p,0 



ob. _ t-3/2 3 p4 -3/2 A ,-2 ,-3/2 ,-3/2 y 



1 



(G) 



The number of radiating electrons is calculated by in- 
tegrating the number density of energetic electrons in- 
side the emitting region, N e (r) — An j r+ArB r 2 n(r)dr ~ 
{(e pl L/Tm p c 3 ){Ar' B /r). Here, Ar/ = Ar' B /T and 
n(r) — Te p in' p (r) are the (observer frame) width and 
number density of radiating electrons inside this region 

3 

By requirement, 7 C > 7 m i n , therefore in calculating 
the observed flux, one can approximate the photon en- 
ergy to be close to The (frequency integrated) 
power emitted by electrons with Lorentz factor 7 m i n is 
-P(7min) = (4g 4 -B' 2 7min)/(9 m e c3 ); therefore the observed 
flux is 

vF ob. = P(^Mr) x p2 

= 2.9 x IO" 7 L 2 rr 3 2 e 2 ^ . 5 eZ£ e B ,_ .5 

x Ar' s , 7 r 2 - 2 dL 2 28 5 3 ,o Co 2 * 2 (P) erg cm- 2 s" 1 , 

(7) 

where — 10 28 5 di, j 28.5 cm is the luminosity distance, 
and a factor P 2 is introduced to transform the result from 
the comoving frame to the observer frame. 

The optical depth is given by r(r) = n' p (r)Ar'crT, 
where the comoving width Ar' and not the comoving 
radiating width Ar^ appears in the equation since elec- 
trons scatter photons outside the radiating region as well. 
The photospheric radius is thus given by 

r(r(r) = 1) = 4^3^7.3 ^ 



r p h 



47rT 3 m p c 3 

= 1.2 x lO^Lr-^Co cm 
The observed synchrotron self absorption energy 
break is calculated usin g standard formula (e.g., 
iRvbicki fc Lightman lll979|) 



31 T 4/5 -8/5 A W o p 



,3/5 r -3/5 



8/5 1/5 



l +2 ^ '13 BJ x 2 

xC /5 Co /5 x(p)eV, 



c e,-0.5 c p/,0 c B,-0.5 



(9) 

where x(p) is a function of the power law index p of 
the accelerated electrons, which is no rmalized to x(p = 
2) = 1. We present in appendix IA. 21 a full derivation of 
this function, and show that its value strongly depends 
on the uncertain value of the power law index p of the 
accelerated electrons. 

While the first four constraints in equation 2] (a-d) are 
common to the majority of bursts, observation of a low 
energy break, which may be attributed to synchrotron 
self absorption frequency is controversial. We therefore 
treat the last constraint in equation ^ separately in sec- 
tion 33 

The constraints set by observations in equation^(a-d) 
can be written with the use of equations ||SJ|-|(HJ in the 
form 

(„) ( l}/ 2 r~ l e 2 f- 2 e 1 ' 2 9 5/2 C 1/2 ^ 2 (v) 

\ a ) \ l + z) ^ '13 e e,-0.5 V,0 e B,-0.5 °p,0 M) * \P! 

= 100ai, 

-2 ^4 -3/2 



<h\ ( 097 ^ T -3/2 3 AJ"' P4, 
y") {T+FJ^ r 13 lAr S,7 i 2 e B,-0.5 ly p,0 



-3/2 >-3/2 
SO 



X^, C0 2 * 2 (P)-1«4, 



100 a 2 , 

1.2Lr2 3 Coa 3> 

e,-0.5 t pL0 tB ^-0.5 U L,28.5 



(c) ri 3 

(d) 2.9L 2 rr 3 2 Ar' B , 7 r 2 - 2 eLn,e;, 1 ne B .-n. s d" 2 



(10) 



3 Due to the requirement r > r p ^, a significant number of pairs 
cannot be created. 
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where the free parameters a\—a^ are introduced in order 
to replace the inequalities in equation ^ by equalities, 
thereby account for the variety of GRB data. 

In order to derive constraints on the values of the free 
model parameters from the set of equations (a-d), 
we note that the parameters e e and es are constrained 
from above by a maximum allowed value of equiparti- 
tion (e e ,-o.5) £s,-0.5 < 1)- The parameter e p i also has 
an upper limit, e p i < 1. Furthermore, the values of 9 p , 
£ and ^(p) (for p > 2) can only be larger or equal to 
unity. In contrast to these constraints, there are no fur- 
ther intrinsic constraints on the values of the isotropic 
equivalent luminosity L, the emission radius r, the fluid 
Lorentz factor T and the comoving decaying length of 
the magnetic field Ar' B . We therefore solve the set of 
equations E3(a-d) to find the values of L, r, T and Ar' B , 
and obtain 

L = 2.7 4,28.5^0.6^,0 Co"' * _1 (P) 

-1/2 1/2 

n 1/1 i \-l j 3/2 _2 1/2 

r 13 = 0.1(1 + 2) dL,28.5e e ;_o.5<v,o e B,-o.5 

K^ ^/2 (p)ar 5 / 4 a l/4 a l/2 ) 

Ar' B , 7 = 4.6 x 10-3 (i + z )-V3 ^ e V» M e -«/a 

T 2 = 3.2 (1 + <28.5 S-^S 4g ^0.5 

(11) 

The values of the free model parameters derived in 
equation 1111 indicate that the prompt emission spectra 
of the majority of the bursts can be explained in the 
framework of the model suggested here. For values of e e 
and cb not far below equipartition and e p i close to unity, 
these results imply that the emission radius should be 
r < 10 12 cm, and that the magnetic field decays on a co- 
moving scale Ar' B ~ 10 4 5 cm. If only w 10% of the elec- 
trons are accelerated in the shock waves, e p i = 0.1, then 
the emission radius is significantly higher, r ~ 10 14 cm, 
and the magnetic field decays after Ar' B ~ 10 6 ' 5 cm. In- 
terestingly, the derived values of the isotropically equiva- 
lent luminosity and the characteristic fluid Lorentz factor 
are not different than their derived values in the standard 
internal shock scenario. We further discuss the implica- 
tions of these results in <21 

3. POSSIBILITY OF A LOW ENERGY BREAK: THE CASE 
OF GRB050820A 

The results obtained in the previous section in equation 
1111 may be applicable to many GRBs that show spectral 
slope vF v oc z^ 4 / 3 below ~ 100 keV. For the majority 
of GRBs observations during the prompt emission phase 
are available only above few keV (BATSE, Beppo-SAX 
or SWIFT BAT-XRT energy range). In most cases, ob- 
servations do not indicate an additional low-energy break 
in the spectrum that might be attributed to synchrotron 
self absorptio n. On the contrary , in some cases (e.g., 
GRB060124: lB,omano et al. 1120061) interpolation of data 
taken in the UV band supports the lack of an additional 
spectral break above ~ 1 eV. 

Even though uncommon to many GRBs, an additional, 
second low energy break may have been obse rved in some 
bursts. In at least one case - GRB050820A ijPaee et al. I 
2005), there are indications for a low energy break at 



< 4keV observed during a gamma-ray/X-ray giant flare 
that occurred 218 seconds after th e burst trigger, and 
lasted 34 seconds ijOsborne 1 12006) . This low energy 
break may be attributed to synchrotron self absorption. 

In order to account for these results in the framework 
of the model presented here, we insert the values of the 
four parameters L, r, Ar' B and T as inferred from the 
observational constraints in eauation llll to the equation 
describing the observed self absorption frequency (equa- 
tion |3- This insertion results in 

F ob. - ^f1 + 7H/5// 5 -29/10 17/5 -7/10 
fc ssa — jb ^t*J U L,28.5 t e,-0.5 t p l,0 fc B,-0.5 

/.-18/5 ~, \ 17/20 -1/4 1/5 1/10 , , 
x0 P X(W"l "2 a 3 a 4 eV ' 

(12) 

where x(p) gives the dependence of £° b sa on the electrons 
power law index p, and is norm alized to xip — 2) — 
1. We present in appendix I A. 21 a full derivation of this 
function, and show there that for values of p in the range 
2 < p < 2.4, this function varies by a factor less than 
4. We can thus conclude that within the framework of 
the model suggested here, the self absorption energy is 
not very sensitive to the uncertain value of the power law 
index p of the accelerated electrons. 

While the the self absorption break calculated in equa- 
tionElis clearly lower than the value of the break energy 
observed in GRB050820A, this equation indicates a very 
strong dependence of the break energy on the uncertain 
values of the parameters e e , e p i and 6 p . The equiparti- 
tion value of e e used in eauation ll2l is an upper limit. If 
the value is e e m 0.1, then the self absorption break is 
observed at ~ 2keV. Similarly, for e p i ~ 0.3 or 9 p w 3 
, the self absorption break occurs at ~ 1 eV. Using the 
results of eauationllllwe find that the values of the four 
other parameters L, r, T and Ar' B are much less sensitive 
to the uncertainties in e e , e p i or 9 p . 

The parameters e e , e p i and 9 p parametrize the post 
shock energy transfer to the electrons, the fraction of the 
electrons population accelerated by the shock waves, and 
the normalized mean random energy gained by proton 
population. All these physical quantities depend on the 
microphysics of energy transfer and particle acceleration 
in shock waves, both of which are not fully understood. 
We cannot therefore, from a theoretical point of view, 
rule out the possibility that the values of e e , e p i and 9 p 
are sensitive to the plasma conditions at the shock form- 
ing region. Equation^] combined with measurement (or 
constraint) on the self absorption frequency, can be used 
to constrain the uncertain values of these parameters. 

4. SUMMARY AND DISCUSSION 

In this work we have presented a model in which the 
magnetic field produced by internal shock waves in GRBs 
decays on a short length scale. Using this assumption, we 
showed that the prompt emission spectra of the majority 
of GRBs can be explained as due to synchrotron radia- 
tion from shock accelerated electrons. We found that 
the required (comoving) decay length of the magnetic 
field is ~ 10 4 5 cm, and that the radiation is produced at 
~ 10 12 cm from the progenitor (eq. 111(1 . These param- 
eter values were found to be relatively sensitive to the 
fraction of electrons population accelerated by the shock 
waves e p i, and can therefore be higher. We showed in JI1 
that the observed synchrotron self absorption energy is 
very sensitive to the uncertain values of the post shock 
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thermal energy fraction carried by the electrons, to the 
mean proton energy 6 P and to the value of e p i, thereby 
argued that the energy of a low energy break is expected 
to vary between different bursts. 

A major result of this work is the characteristic de- 
cay length of the magnetic field deduced from observa- 
tions, ~ 10 4 5 cm. This value is significantly shorter 
than the standard assumption used so far, that the mag- 
netic field strength is approximately constant through- 
out the comoving plasma width, w 10 10 — 10 11 cm. 
Still, the electron crossing time of the magnetized re- 
gion is long enough to allow electrons acceleration to 
high energies. Equating the electron acceleration time, 
t acc ~ 7m e c 2 /(cgB') and the electron crossing time, 
Ar' B / c, gives an upper limit on the electron Lorentz fac- 
tor, 7maXil = (Ar' B qB')/(m e c 2 ) ~ 2 x 10 7 Ar' B A 5 B' 6 , 
where B' — 10 6 B' 6 G. This value is larger than the 
maximum electron Lorentz factor obtained by equat- 
ing the acceleration time and the synchrotron cooling 

time, 7max , 2 = (3/2)(m e c 2 )/(g 3 J B') 1/2 ~ 10 5 ^" V2 . We 
thus conclude that within the magnetized region, for 
e B ^ 10~ 3 an upper limit on the accelerated electron 
energy is set by the synchrotron cooling time, and not 
by the physical size of this region. 

The parameters a\ — introduced in equations I1UI 
^] account for the difference between the variety of GRB 
data and the characteristic values considered in the an- 
alytical analysis. As a concrete example, cooling energy 
e° b ' larger than 100 keV is accounted for by considering 
Q!2 > 1, which implies through equation 1111 that high 
isotropic equivalent luminosity is required in order to ac- 
count for an observed flux ~ 10 _7 ergs _1 cm -2 . This re- 
sult is understood as due to the low radiative efficiency 
in the case e° b - > 

The ratio found here between the decay length of the 
magnetic field and the comoving shell thickness, 

§^ = 1.5xlO- 5 e-i . 5 e p; ,oe^-o. 5 ^o*- 1 (p) 

1/2 -1/2 

xa-,/ a 2 ck3, 

(13) 

is based on fitting the GRB prompt emission spectra. 
In earlier work, bas ed on modeling afterglow emission 
ijRossi fc Rees ll2003j) . this value was thought to be too 
low, due to the high ambient medium density it implies 
during t he afterglow emission phase. However, in the 
work by IRossi &; Rees I l|2003|) detailed modeling of af- 
terglow data was not performed, due to lack of avail- 
able data. Moreover, the well established connection be- 
tween long GRBs and core collapse of massive stars (e.g., 
iPe'er Sz Wiiers I |2006t and references therein) indicates 
that indeed the ambient medium density may be higher 
than previously thought. 

The values of the parameters found in eauation llll im- 
ply that the comoving number density of the shocked 
plasma (equation |2J is 

n' p = 5.0 x 10" (1 + z) 4 /3 d£% e ;% «{g e B %, 
x9;l *- 7 /3 (p ) a 3 / 2 a- 1/e a- 2/3 a^ 1 ' 3 cm" 3 . 

(14) 

For this value of the comoving number density, the 



plasma skin depth is given by 

\ ~ _ n o f-i , -A-2/3 jl/3 5/3 -11/6 1/3 

A - ~ V~*{L + Z) d L,28.5 e e,-0.5 e pl,0 e B-0.5 

x 6 2 p * 5 / 3 (p) a- 3/A al /12 a\ /3 a\ /e cm, 

(15) 

where uj pe = (A-Kq 2 n' p /m e ) 1 / 2 is the plasma frequency. 
This value of the skin depth implies that the magnetic 
field decays on a characteristic length scale 

= 2xl0^(l + z)-/3^, 5 ^ 5 «^, 
*e-l^{ P )a-^a-^ 2 ay 3 aT 

(16) 

skin depths. This decay length of the magnetic field is 
four orders of magnitude shorter than the characte ristic 
scale « 10 9 skin depth assumed so far l)Piran 1120051) . On 
the other hand, it is three orders of magnitude longer 
than the maximum length scale of magnetic field genera- 
tion that can be calculated usi ng state of the art n umer- 
ical models llSilva et al. 1120031 IFrederiksen et al. 1120041 
iNishikawa et al. Il2005|) . The results obtained here are 
based on the interpretation of GRB prompt emission 
spectra. They can therefore serve as a guideline for the 
characteristic scale needed in future numerical models of 
magnetic field generation. 

The results presented in equations ^2 113 an d GUI in- 
dicate that the value of e e should be close to equipar- 
tition. The value of cb on the other hand, is less con- 
strained, and values as low as 1-2 orders of magnitude be- 
low equipartition are consistent with the data (the strin- 
gent constraint on the value of es is obtained by the self 
absorption energy, equation I12fl . The results presented 
in eouation ll II indicates that a low value of e p i results in 
large emission radius and large decay length of the mag- 
netic field. Thus, low value of e p i implies that the model 
presented here can account for late time flaring activities 
observed in many GRBs, that may originate from shell 
collisions at large radii. A lower limit on the value of e p i 
can be set by the requirement that the emission radius 
is not larger than the transition radius to the self similar 
expansion, ~ 10 16 cm which marks the beginning of the 
afterglow emission phase. From this requirement, one 
obtains e p i > 10~ 2 . 

Generation of magnetic field and particle ac- 
celeration in s hock waves are most probab ly re- 
lated issues dKazimura et a lA Il998t iSilva et al.\ 
20031: IFrederiksen et al. I 120041: Ififededal et al. I 120041 
Nishikaw a et al. 1 12005'). We therefore anticipate that 
the answers to the theoretical questions raised by the 
model presented here, about the requirement for high 
values of e e and es, the uncertainty in the value of e p i 
and the characteristic decay length of the magnetic field, 
are related to each other. 

An underlying assumption in the calculations is that 
the values of the free parameters are (approximately) 
constant inside the emitting region. In reality, this of 
course may not be the case. We introduced here a new 
length scale Ar' B , characterizing a length scale for the 
decay of es- It can be argued that within the context of 
this model the decay length of e e is not shorter than Ar' B . 
However, by defining Ar' B as the shortest length within 
which both e e and eb maintain approximately constant 
values, the results presented here hold. 

The emission radius r « 10 12 cm found, implies that 



(> 
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this model can account for observed variability as short 
as r/r 2 c « 1 ms. The observed GRB prompt emission 
spectra are usually integrated over a much longer time 
scale, of few seconds. This can be accounted for in our 
model, cither by assuming low value of e p i, or by adopt- 
ing the commonly used assumption that the long dura- 
tion emission is due to extended central engine activity, 
which continuously produces new shock waves and re- 
freshes existing shock waves. 

The results presented here are applicable to a large 
number of astrophysical objects, in which magnetic field 
generation and particle acceleration in shock waves are 
believed to play a major role. Such is the case for the 
study of afterglow emission from GRBs as w ell as emis- 
sion from supernovae remnants (see, e.g., iChevalierl 
1992, for the case of SN1987A). Additional astrophys- 
ical sources in which strong shock waves and magnetic 



fields occur are a ctive galactic nuclei (AGNs) and jets 
in micro-quasars ijFender [ [20031 . Current observational 
status of these objects confines sy nchrotron emitting re - 
gions only on a scale of ~ 10 13 cm l|Dhawan et al. 1 2000'). 
If the length scale of the magnetic field inferred from ob- 
servations in these objects is found in the future to be 
similar to the value found here, i.e., « 10 5 skin depths, 
this may serve as a strong hint toward understanding 
magnetic field generation in shock waves. 
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APPENDIX 

THE DEPENDENCE OF THE BREAK ENERGIES ON THE POWER LAW INDEX P OF THE 

ACCELERATED ELECTRONS 

Electrons minimum Lorentz factor, "f m in 

We assume that a fraction e p i of the electrons are accelerated to a power law energy distribution p above 7 m ; n and 
below 7 m ax- While the initial value of 7 m i n depends on the bulk Lorentz factor of the flow, a quasi steady state of the 
electron distribution is formed, in which the value of 7 m i n depends on the number and energy densities of the acclerated 
particles. As this happens, the Lorentz factor 7 m ; n can be calculated given the number and energy densities of the 
accelerated electrons. The electron energy distribution is given by dn/dj = Aj~ p , where A is a numerical constant. 
Integrating this function relates the values of 7 m i n and A to the number and energy densities of the energetic electron 

component, e pi n el ~ e pi n' p = J^(dn/d-/)d-y = A(l - p)- 1 - 7^7n) > and 



dn 



-jdj 



rripC 2 x 



Dividing u e i by e p in e im e c 2 eliminates A from the equations, 



A log 

A 

(2-P) 



7m 
I max 



7 2 " p 

/ mm 



(p = 2), 
(p * 2). 



(Al) 



(P = 2), 



We can now write the value of 7 
^(p) is given by 



in in 7min 
1 



*(P) 



(A2) 

u e i/(m e c 2 e pl n e i)]{l - 7min/7max) log(7 ma x/7min) _ " 1 x *(p), where 

(A3) 



(P = 2) 
&> ^ 2). 



The function ^(p) is plotted in figure |A"T1 for two representative values of (7min/7max) 4 - 

Self absorption energy, e°^ a 

The synchrotron self absorption coefficient for a power law dis tribution of electrons with p ower law index p radiating 



ine synenrotron sell absorption cocmcicnt lor a power law distribution 01 electrons witn pc 
in magnetic field B' is calculated using standard formula fe.g.. lRvbicki fc Light man Ill979|) . 



27 9 4 



167r 2 mic 5 V 27rm 3 c 5 



3q 



p/2-1 



r(f)r(I) 



- ) A R'P/2+l l/ -(p/2+2) 

3, 



where the constant A is calculated using eauation lAll 



A = u el x 



lo 



(P ^ 2), 



'"2)7^ 



(A4) 



(A5) 



4 The value of 7 max can in principle be found from physical constraints on the acceleration time. However, we find this way of presentation 
much clearer. 




which, upon insertion of 7min can be written as 

-l 



A = u ei log 



In 



7d 



U e l 



e p in e im e c^ 



p-2 



where 



1 (Tmin/Tmax) 

(P = 2), 



v^- 2 (p)£(p), 



(A6) 



(A7) 



Inserting the numerical values of the magnetic field and the peak frequency v pea k — e^ /Th (see equations |21|5jl into 
the self absorption coefficient equation IA4I using the value of A found in equation IA6I one obtains the synchrotron 
self absorption coefficient at the peak frequency, 



1.00 X 1U L r 13 1 2 e e , — 0.5 pl,Q B, — 0.5 p.O '■O 

2-1 >(¥M^)" 



x(3.0 x 10 12 ) p / 



^(p)$- 6 (p) cm" 1 . 



(A8) 



r(l)r(I) 

The self absorption optical depth t v = Ar' B a v is smaller than unity at v = v pea k- Since (by demand) the electrons are 
in the slow cooling regime (i.e., 7 m ; n < 7 C ), the power radiated per unit energy below e m — /V is proportional to 
(e/em) 1 / 3 , and the energy below which the optical depth becomes greater than unity, e ssa = s m T^= Up(;ak , is 

-4/5 ^4/5 

(A9) 

mm 



-ob. _ _31_ t 4/5 -8/5 
-ssa — l+ 2 u ' 13 



,3/5 p-3/5 
B,7 1 2 



X(3.0X l 12)(3/5)( P /2-l) 



-1 8/5 1/5 a-*/" r- 
e,-0.5 fc p;,0 fc B,-0.5 ^p.O ^0 

3p + 2 -jp^ 3P + 22 N 1 3 / 5 
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(compare with equation [5J . Therefore, the definition of the function x(p) is 



X {p) = (3.0 x io 12 )(3/5)(p/ 2 -D 

This function is plotted in figure ^2] 

Inserting the parametric dependence on the value of of the four parameters found in equation into equation 
IA9I leads to x(p) — x(p) x ^~ 19 / 10 (p)- Graph of this function appears in figure 1X51 Note that while x(p) shows a 
very strong dependence on the value of p, the function x(p) varies by a factor less than 4 in the range 2 < p < 2.4. 
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